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Warming in complex physical systems, in particular global warming, attracts significant con- 
temporary interest. It is essential, therefore, to understand basic physical mechanisms leading to 
overheating. It is well known that application of an electric field to conductors heats electric charge 
carriers. Often an elevated electron temperature describes the result of the heating. This paper 
demonstrates that an electric field applied to a conductor with discrete electron spectrum produces 
a non-equilibrium electron distribution, which cannot be described by temperature. Such electron 
distribution changes dramatically the conductivity of highly mobile two dimensional electrons in a 
magnetic field, forcing them into a state with a zero differential resistance. Most importantly the 
results demonstrate that, in general, the effective overheating in the systems with discrete spectrum 
is significantly stronger than the one in systems with continuous and homogeneous distribution of 
the energy levels at the same input power. 



The energy distribution of hot electrons in response 
to an electric field is a long-standing, interesting prob- 
lem in condensed matter and plasma physics !lj. Despite 
the simplicity of the experiments, interpretation of re- 
sults is based significantly on theoretical suggestions. It 
is widely accepted that overheated, degenerate electron 
systems are well described by an elevated electron tem- 
perature Te- This approach is based on an analysis of the 
Fokker-Plank type equation with almost elastic scatter- 
ing y. Most dramatically the temperature broadening 
of the electron distribution function affects kinetic prop- 
erties of electron systems with a discrete spectrum, in 
particular, electrons in a magnetic field. Shown by L.D. 
Landau [3| a quantization of the electron motion in a 
magnetic field results in the discrete spectrum of two di- 
mensional electrons. Direct consequences of the Landau 
quantization are quantum magnetic oscillations in met- 
als: Shubnikov-de Haas effect and de Haas-van Alphen 
effect y, |5[. At the thermal equilibrium the amplitude 
of the quantum oscillations depends exponentially on the 
electron temperature Tg. Often the oscillations are used 
as an electron thermometer. 

Below we show that temperature Te is not an appro- 
priate parameter for 2D electrons overheated by a direct 
current (DC) in magnetic field. In other words the non- 
equilibrium distribution function of the DC biased 2D 
electrons is not described by a temperature. The DC 
bias originates a non-uniform spectral diffusion of the 
2D electrons over the discrete Landau spectrum result- 
ing in a non-equilibrium electron distribution oscillating 
with the energy. Our finding is in agreement with recent 
theoretical consideration of the effect of the DC electric 
field on the conductivity of 2D electrons [6| and recent 
experiments in this field [7J, [sj, |9J, llO| . 

Our samples are cleaved from a wafer of a high- 
mobility GaAs quantum well grown by molecular beam 



epitaxy on semi-insulating (001) GaAs substrates. The 
width of the GaAs quantum well is 13 nm. Three sam- 
ples (N1,N2,N3) are studied with electron density ni=8.2 
xlO^^ (m-2), 712=8.4x10^5 (ni-2), n3=8.1xl0^'^ (m^^) 
and mobility /ii=85 (m^/Vs), /i2=70 (m^/Vs) and ^J,3=82 
(m^/Vs) at T=2K. Measurements are carried out be- 
tween T=0.3K and T=10 K in magnetic field up to 1.5 
T on d=50 iim wide Hall bars with a distance of 250 fj,m 
between potential contacts. The differential longitudinal 
resistance is measured at a frequency of 77 Hz in the lin- 
ear regime. Direct electric current (dc bias) was applied 
simultaneously with ac excitation through the same cur- 
rent leads (see insert to Fig. 1). All samples demonstrate 
similar behavior. Below we present data for sample Nl. 
Figure la shows the longitudinal resistance r^x = 
dVxx/dl [ll| as a function of the magnetic field B at 
zero DC bias and at different temperatures as labeled. 
At B >0.4 (T) the resistance exhibits the Shubnikov 
de Haas oscillations. The amplitude of the oscillations 
decreases significantly with temperature increase. The 
temperature reduction of the SdH amplitude is result of 
an effective (exponential) cancellation of the periodic os- 
cillations of the electron conductivity a{e) with energy 
e weighted with the first derivative of the equilibrium 
(Fermi-Dirac) distribution function df /de and averaged 
over all possible energies [5|: 



a{e){-df/de)de 



(1) 



Figure lb shows the longitudinal resistance rxx = 
dVxx/dl as a function of magnetic field B taken at three 
different dc biases as labeled at a lattice temperature 
T=2.1K. The dc electric field changes significantly the 
longitudinal resistance in magnetic fields above 0.15(T). 
At i? > 0.15 (T) the disorder broadening Ti/rq of the 
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FIG. 1: (Color online) (a) Dependence of longitudinal resis- 
tance r^x on magnetic field B with no dc bias Idc applied at 
different temperatures as labeled. Insert presents the experi- 
mental setup, (b) Dependence of longitudinal resistance r^x 
on magnetic field B at different dc bias Idc as labeled. Insert 
presents the bias dependence at fixed magnetic field 5=0.797 
(T). Temperature is 2.1 (K). 



Landau levels is less than the Landau level separation 
tiujc and the electron energy spectrum (density of states) 
becomes modulated. The modulation increases with the 
magnetic field. At S ^ 1 (T) the Landau levels are 
completely separated from each other and the electron 
spectrum is discrete. 

The insert to figHh demonstrates a typical dependence 
of the longitudinal resistance r^x on the dc bias in a mag- 
netic field {B — 0.797 (T)), corresponding to an oscilla- 
tion maximum. For small dc bias the resistance decreases 
approaching zero. As the dc bias is raised further, the 
resistance exhibits a reproducible sharp negative spike at 
/tft=8.4 (/iA) and, then, stabilizes near zero at higher bi- 
ases. The zero differential resistance state is observed in 
a broad range of magnetic fields at temperatures below 
few Kelvins[9|, llO|. 

The curves presented in figUJa reveal a dramatic and, 
in many respects, unexpected evolution of the SdH os- 
cillations with the dc bias. The behavior is significantly 
different from the temperature evolution shown in fig. la. 
Both oscillating and smooth parts of the dependences 



FIG. 2: (Color online). Longitudinal resistance Vxx plotted 
vs inverse magnetic field 1/B at different dc bias as labeled. 
Temperature T=2.1 (K). 



in fig. lb exhibit substantial changes. The smooth part 
of the resistance displays a strong reduction with the dc 
bias in magnetic fields B > 0.15 (T), at which the density 
of electron states oscillates with energy. The oscillating 
part of the resistance reveals a rather complicated behav- 
ior with the dc bias. 

Figure 2 presents the bias induced evolution of the 
quantum oscillations in greater detail. Plotted as a func- 
tion of the inverse magnetic field 1/B the oscillations are 
periodic [J| . Although the fundamental periodicity is the 
same for all biases, the form of the oscillations and/or 
the harmonic content are changed significantly with the 
bias. The most dramatic changes occur at maxima of the 
oscillations. With the increasing bias the maxima evolve 
gradually into minima with essentially zero differential 
resistance. On the other hand the minima tend to form 
maxima and, finally, separate the zero resistance regions 
at the end of the evolution. For further discussion it is 
important to realize that the high frequency content of 
the oscillations with respect to the fundamental period- 
icity increases with the dc bias increase. 

Figure 3a demonstrates quantitatively the evolution of 
the first and second harmonics of the quantum oscilla- 
tions with temperature at zero bias. At thermal equilib- 
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FIG. 3: (Color online), (a) Fourier spectrum of the quan- 
tum oscillations at different temperature as labeled. Insert 
shows temperature dependence of amplitude of first (at 17.1 
T) and second (at 34.2 T) harmonics of the quantum oscil- 
lations. Zero dc bias, (b) Fourier spectrum of the quantum 
oscillations at different dc bias as labeled. Insert shows the 
bias dependence of amplitude of first (at 17.1 T) and second 
(at 34.2 T) harmonics of the quantum oscillations. Lattice 
temperature is 2. IK. 



rium the distribution function / does not oscillate with 
energy and, therefore, the harmonic content of the quan- 
tum oscillations is determined only by the harmonic con- 
tent of the conductivity a{e) oscillating with energy e 
due to Landau quantization. The amplitudes of the first 
A^(T) and the second A2ui{T) harmonics exhibit expo- 
nential decrease with the temperature due to the spectral 
average of the oscillations of the conductivity a(e) over 
the range of energies Ae corresponding to the thermal 
broadening of the distribution function Ae ~ kT. At 
higher frequency the spectral cancellation is more effec- 
tive and, as a result, the higher harmonics of the quantum 
oscillations {A2uj{T)) are exponentially smaller than the 
amplitude of the fundamental periodicity Ai^{T). The 
insert to figinh- demonstrates that the amplitudes of the 
first A^ and the second A2uj harmonics as well as their 
ratio R = A2^/A^ decrease significantly with the tem- 
perature and approach the zero value asymptotically. 

Figure 3b demonstrates the effect of the dc bias on the 



spectrum of the quantum oscillations. Notably differ- 
ent behavior of the first and, especially, the second har- 
monics with dc bias is observed. The first (fundamental) 
harmonic does not approach the zero value asymptoti- 
cally. Instead, as it is shown on the insert to fig. 3b, the 
amplitude crosses the zero value and becomes negative at 
Idc = 7.5 fiA. In other words, the phase of the fundamen- 
tal oscillations changes by 180". This corresponds to the 
evolution presented in figure 2b, where maxima become 
minima with the bias increase. Even more unusual is the 
behavior of the second harmonic of the oscillations. One 
can see that the magnitude of second harmonics does not 
exponentially decrease with the dc bias contrary to the 
usual expectation. Instead the second harmonic crosses 
the zero and reaches a value, which is even bigger than 
the amplitude ^21^(0) at zero bias. At the same time the 
ratio R between the amplitudes of the second and first 
harmonics is on the order of unity, which, again, cannot 
be explained by the temperature broadening of the distri- 
bution function. Thus the results indicate that the effect 
of the direct current on resistance of the 2D electron sys- 
tem with discrete energy spectrum cannot be described 
by an increase of the electron temperature. 

The enhancement of the higher harmonics of the quan- 
tum oscillations with the DC bias follows directly from 
eq.(l) assuming additional, bias induced oscillations of 
the non-equilibrium distribution function with the en- 
ergy. Such oscillations have been found theoretically in 
DC biased 2D electron systems recently^]. Due to the 
conservation of total electron energy in the presence of 
an electric field E and elastic electron-impurity scatter- 
ing the kinetic energy of an electron depends on electron 
position f: e{r) = eo + eEr. As a result, diffusion motion 
of the electron in real space originates a diffusion of the 
electron in energy space. The coefficient of the spectral 
diffusion Z?e is proportional to the coefficient of the spa- 
tial diffusion D: D, = {cEfD - (Sr)"^ . The spatial (and 
so spectral) diffusion is most effective in the center of the 
Landau levels, where the density of states is high, grad- 
ually decreases away from the center and is suppressed 
considerably between Landau levels, where the density of 
states is small. 

The spectral diffusion generates an electron spectral 
flow Jf from low energy regions (occupied levels) to high 
energies (empty levels). The spectral flow is proportional 
to Df^ and to the gradient of the distribution function 
df/de: J^ = D{e)-df/de. In a stationary state the spectral 
electron flow J^ is constant. As a result, the gradient of 
the distribution function df/de is strong in the regions 
of weak spectral diffusion (between Landau levels) and is 
small in the regions with strong spectral diffusion (centers 
of the Landau levels). 

A weak inelastic relaxation of the distribution func- 
tion with an inelastic scattering rate l/ri„ exists at low 
temperature [,6|. The weak inelastic scattering does not 
change considerably the spectral diffusion. However on 




and/or density of states makes the net conductivity cr 
(see eq.l) significantly smaller than the unbiased value. 
This explains the strong reduction of the conductivity 
(and resistivity r^x III) with the DC bias at magnetic 
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FIG. 4: (Color online), (a) right axes: density of states of 
2D electrons in magnetic field B=0.797 (T) vs energy; left 
axes: numerical calculations of electron distribution function 
modified by the spectral diffusion. The following parameters 
are used in the simulation: magnetic field B=0.797 T, tem- 
perature T=8 K, dc electric field E — 1003 (V/m), disorder 
broadening of Landau levels h/Tq=h/1.9 (1/ps) and inelastic 
scattering time Ti„ — 150 (ps). The times are obtained by 
methods presented in paper [7] ; (b) gradient df/de of the dis- 
tribution function modified by the spectral diffusion shown 
above and the gradient of the thermal equilibrium function 
corresponding to different temperatures as labeled. 



long times t ^ Tin the inelastic processes stabilize the 
bias induced deviations of the distribution function from 
the thermal equilibrium. It is taken into account below 
in a numerical simulation of the spectral diffusion. The 
spectral diffusion equation (see eq. (9) in [6,]) is solved 
numerically. 

Results of the numerical simulation of the spectral dif- 
fusion are shown in figure 4a. As expected the spectral 
diffusion produces periodic variations of the distribution 
function /(e) with energy. It is clear, that an elevated 
temperature Tg cannot describe this oscillating behavior. 
The additional oscillations of the distribution function 
enrich significantly the harmonic content of the prod- 
uct a{e) ■ {—df/de) and, therefore, in accordance with 
eq.(l), the spectrum of the quantum oscillations. Fur- 
thermore shown in fig. 4b the considerable decrease of the 
value df/de at the maximums of the conductivity cr(e) 



field B > 0.15 T shown in fig. lb. 

The strong reduction of the gradient of the electron 
distribution function df(e)/de near the maximums of the 
conductivity a{e) and/or the density of states imitates a 
strong overheating in these energy regions. Indeed fig. 4b 
demonstrates that at a maximum of the density of states 
located at the Fermi level /ii e — /i = 0, the gradient of the 
DC biased distribution function df{e)/de equals the one 
corresponding to the thermal equilibrium at temperature 
of 32 K. Thus at the conductivity maximum a{e = fi) the 
effective overheating of the 2D electrons (32K) is signif- 
icantly stronger the actual additional broadening of the 
distribution function (^ 1 — 2K) induced by the DC bias. 
This provides the dramatic change of the transport prop- 
erties of the system (see fig. lb). 

The strong reduction of the df/de is result not only 
of the high spectral diffusion near the spectral maxima, 
but, mostly, is a result of the weak spectral diffusion 
between them. In fact the weak diffusion between the 
spectral maxima controls the global spectral flow J^. In 
some sense the regions of the phase space with the weak 
spectral diffusion provide a "thermal" isolation of the re- 
gions with the strong spectral diffusion. This DC bias 
induced, thermal stratification of the phase space facili- 
tates the strong overheating in the regions essential for 
the net electron transport. Similar spectral diffusion and 
the thermal stratification are expected in 2D electron sys- 
tems affected by an electromagnetic radiation. In this 
case the spectral diffusion is a result of quantum transi- 
tions between energy levels induced by the external radia- 
tion. This points toward a generic nature of the phenom- 
ena. Thus the presented results suggest, that in systems 
with a discrete or strongly modulated energy spectrum 
the effective overheating is significantly stronger than 
the one in systems with weakly modulated, homogeneous 
spectrum at the same input power. 
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